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ABSTRACT 

We present L' band (3.8 nm) MMT/Clio high-contrast imaging data for the 
nearby star GJ 758, which was recently reported by Thalmann et al. (2009) to 
have one - possibly two- faint comoving companions (G J 758B and "C" , respec- 
tively). GJ 758B is detected in two distinct datasets. Additionally, we report a 
possible detection of the object identified by Thalmann et al as "GJ 758C" in 
our more sensitive dataset, though it is likely a residual speckle. However, if it is 
the same object as that reported by Thalmann et al. it cannot be a companion 
in a bound orbit. GJ 758B has a H-L' color redder than nearly all known L-T8 
dwarfs. Based on comparisons with the COND evolutionary models, GJ 758B 
has Tg ~ 560 K-fo^^ and a mass ranging from ~ 10-20 Mj if it is ~ 1 Gyr old to 
~ 25-40 Mj if it is 8.7 Gyr old. GJ 758B is likely in a highly eccentric orbit, e ~ 
0.73-0-21, with a semimajor axis of ~ 44 AU-'^^au _ Though GJ 758B is sometimes 
discussed within the context of exoplanet direct imaging, its mass is likely greater 
than the deuterium-burning limit and its formation may resemble that of binary 
stars rather than that of jovian-mass planets. 

Subject headings: planetary systems, stars: brown dwarfs, techniques: high an- 
gular resolution 



Introduction 



High-contrast imaging surveys have recently identified many faint, cool companions to 
nearby stars whose inferred masses are between that of Jupiter and the deuterium burning 
limit (~ 13 Mj). The companions to nearby A stars - HR 8799, Fomalhaut, and {3 Pic 
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- orbit at separations less than ~ 100 AU, have small ma ss ratios, and genera . 
scale d -up versions of gas gia nt planets in our solar system (IMarois et al.ll2008a ;_ 



2008a: iHinz et al 



Metchev et al. 



20l31currie et al 



20081: iLagrange et al.M20101). with tempera tures that are (likely) < 1500 K (IMarois et al- 



ly resemble 



Kalas et al. 



20 



Oal). comparable to that for many L and T dwarfs 



20061 : Iheggett et al.lboiol ). 0. Other directly imaged companions have 



(e.g. 

similar temperatures but orbit low- mass stars and brown dwarfs, typically at wider sep 



arations and/or with much larger mass ratios fe.g. IChauvin et al.l |2004J : lltoh et al.l 12005 



Luhman et al.l 120061 : iLafreniere et al.l l2010l : iTodorov et al.l 120101 ) , indicating that they plau- 
sibly represent the low-n aass tail of objects formed by molecular cloud core or protostellar 



disk fragmentation (e.g. iLodato et al.l 120051 ). The planetary companions to A stars pose 



strong challenges for even the most efficient core accretion models of gas giant plan et for- 



mation, given the difficulty of forming massive cores at 10-100 AU (e.g. lRafikovll2010l ) prior 
to the dispersal of the protoplanetary disk. The configuration of at least one system with a 
low-mass primary and planet-mass companion, 2MASS J04414489+2301513, implies that the 
fragmentation of molecular cloud cores can prod uce objects < 5-10 M j, below the classical 
"opacity-limited" minimum fragmentation mass ( iTodorov et al.ll2010l ). 



The faint companion and candidate companion to the nearby star GJ 758 reported by 
Thalmann et al. (2009; hereafter T09), GJ 758B and "GJ 758C", present an intriguing 
contrast to planets orbiting A stars and to other substellar-mass L/T dwar fs. GJ 758 is 



much later in spectral type (G8V) and lower in mass (M^ = 0.97 ± 0.03 M©, iTakeda et al. 
2003) than HR 8799, Fomalhaut, and (3 Pic (1.5-2.1 M©). On the other hand, the projected 
separation for GJ 758B (1.858", PA = 198.18° on 2009 Aug. 6, or ~ 29 AU) and "GJ 758C" 
(1.118", PA = 219.1 6°, or ~ 18 AU) are comparable to the se parations for planets orbiting 



HR 8799 an d 3 Pic (IMarois et al.ll2008al : iLagrange et al.ll2010l ) and within the brown dwarf 



desert (e.g. iKraus et al.l l2008l ) . Based on GJ 758B's H band magnitude and age (0.7-8.7 
Gyr), T09 argue that the companion likely has a temperature of 549-637 K (M = 11.7- 
48.5 Mj), making it the coldest known companion to a Sun-like star. Since T09 image 
the system only in H band, the companion's temperature -and thus luminosity - remain 
observationally unconstrained. Inferred masses critically depend on these properties, so new 
longer wavelength data help constrain whether GJ 758B better resembles the low mass-ratio 
planetary companions to A stars or high mass-ratio binary star companions to lower-mass 
stars. 

Here we report the L' band detection of GJ 758B and candidate detection of "GJ 758C" 
with the Clio camera at the 6.5m MMT telescope. Combined with the H band photometry 



^Dynamical constraints for the masses of some of these companions - e.g. those orbiting HR 87 99 - are 
consistent with luminosity-derived estimates ( Chiang et al. 2009 : Fabrvckv and Murrav-Clav 2010[) 
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from T09, we estimate the temperature, luminosity, and mass of GJ 758B. Combining all 
astrometric datapoints for GJ 758B, we determine its range of allowed orbits. 



Observations and Data Reduction 



GJ 758 was imaged under photometric conditions on May 27, 2010 and June 2 , 2010 



at the 6.5m MMT tele scope with the upgraded Clio mid-IR camera (IHinz et al.l 12006 



Sivanandam et al. 2006). As described in online documentatioij^, the new Clio detector 
(Clio-2) operates from 1.65 fim to 4.8 /im, covering standard Mauna Kea filters and narrow- 
band filters centered on wavelengths between 3 and 4 /im. Both sets of data reported here 
were o btained in the L' fi lter (3.8 fim). All data was taken in angular differential imaging 
mode (IMarois et al.l 120061 ). which keeps the instrument rotator fixed, allowing the field of 
view to rotate with time. The May data consist of coadded frames of 9s each for a total 
integration time of 2520s. The data were taken through transit, yielding a total field rotation 
of 192 degrees. The June data consist of coadded frames of 9.6s each for a total integration 
time of 1536s but were obtained well after transit for a total field rotation of 16.7 degrees. 
In both datasets, the star was dithered along the detector by 5-10" every 48-90s for sky 
subtraction. For precise astrometric calibration during our run, we observed the double stars 
HIP 88817/88818 and HD 223718. The detector orientation is offset by 2.53° ± 0.15° coun- 
terclockwise from true north along the y axis; the pixel scale is 0.029915 ± 8xl0~^ "/pixel. 
Independent astrometry performed by multiple coauthors confirm these values. Based on 
this pixel scale, our image FWHM is ~ 0.15" for both datasets. 

Our image processing method closely followed the ADI/LOCI (Angular Diffe rential 



Imaging/Loc ally- Optimized Combination of Images) reduction procedure described by lMarois et al, 
(120061, |2008aJ) and iLafreniere et al.l (120071 ) used to maximize the companion signal relative 
to speckle-dominated noise. The subtracted residual images produced by our ADI/LOCI 
pipeline were then derotated, median combined, and convolved with a gaussian kernel equal 
to the image FWHM to produce final science images. Separately from our LOCI reduc- 
tion, we performed a simple ADI reductio n by derotating, high-pass filtering, and median 
combining the PSF subtracted images (e.g. iHinz et al.ll2010l ). 



Figure |T] shows the final reduced images for the May dataset (left) and June dataset 
(right) produced with the LOCI pipeline. The May dataset clearly reveals a point source 
at a separation of 1.823" ± 0.015" (a 



28.4 AU) and position angle of 199.76 ± 0.15 
degrees. Within astrometric errors, we identify a point source at the same position in the 



proj 



^http://zero. as. arizona.edu/wiki/doku.php 
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lower-quaility June data (1.827" ± 0.043"|j), indicating that its detection is secure. We 
recover the point source in our more simple ADI reductions as well. This separation is 
comparable to that for GJ 758B as reported by T09 (1.858" ± 0.005"); the difference in 
position angle (5(PA) = 1.58") is consistent with an object undergoing counterclockwise 
motion with respect to the primary. The velocity in the plane of the sky of GJ 758B, derived 
by comparing our May 2010 position with previously reported positions (T09), is 1.3 ± 0.2 
AU/yr. This is less than the minimum escape velocity at this projected separation (1.64 
± 0.03 AU/yrx(M^/O.97M0)i/2x(ap/28.5O AUy^^'^), suggesting a bound orbit. Thus, we 
conclude that we have detected GJ 758B. 

The LOCI reduction of the May data also reveals a second, candidate point source at 
1.101" ± 0.015", comparable to that reported for the candidate companion, "GJ 758C", in 
T09. However, the large difference in position angle (~ 20°) would imply a space velocity of 
7.9 ± 0.3 AU/yr, nearly four times greater than the escape velocity of 2.08 AU/yr. Thus, if 
our candidate point source is "GJ 758C", it cannot be in a bound orbit around the primary. 

To assess the significance of the GJ 758B detection and candidate "GJ 758C" detection, 
we compute the standard deviation and signal-to-noise of pixel values in concentric annuli 
(e.g. T09). The signal-to-noise is ~ 6.4 at GJ 758B's position and ~ 3 at the candidate 
"GJ 758C" position. Thus, the signal-to-noise of our possible "GJ 758C" detection is only 
marginally significant (~ 3cr), even though it is locally well separated from large background 
fluctuations. Bona- fide detections in speckle- noise limited regions typic ally need to be greate r 



than 5a in order to rule out false "detections" from residual speckles (iMarois et al.ll2008bl ). 
Because of the chance that "GJ 758C" could be a residual speckle and the lack of a detection 
in our June data, we consider the second point source in our May data to only be a candidate 
detection of the point source identified by T09 as "GJ 758C" □ Since only the detection of 
GJ 758B is secure, we focus on it in our analysis. 



■^The much larger error bars result because the companion is contaminated by some residual speckle noise 

''Chance superpositions of substellar objects close to stars like GJ 758 with L' band brightnesses similar 
to our candidate detection are unlikely t o be frequent. As s umin g a flat substellar mass function (a — 0) with 
a space density described in Table 4 of iBurgasser et al.l (120041) using the abs olute magnitude vs. spectral 
type of L/T dwarfs at the Spitzer/IRAC [3.6] band from lLeggett et aL ( 2010l ) as a proxy for L' brightness, 
the probability of contamination by a background brown dwarf with Te =500-1500 K is < lO^'^ %. 
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3. Photometric and Astrometric Analysis 



Photometry for GJ 758B from the May data was performed with IDLPHOT, using a 2.5 
pixel aperture radius and a background annulus between 2.5 and 5 pixels. In all exposures, 
the stellar PSF core is saturated. For photometric calibration, we compare the GJ 758B flux 
to that for HD 223718A , which was observe d immediately after GJ 758. HP 223718A is a n 
F5V star (Te ~ 6440 K. ICurrie etaDboiObh and should have K-L' = 0.04 jTokunagalboOof ). 
Based on its 2MASS photometry (Ks= 6.51) and using t he updated versi on of the 2MASS 
color transformations between 2MASS and CIT systems (jCarpenteii 1200 ll ) . its estimated L' 
magnitude is 6.49. Both the ADl and LOCI reduction procedures attenuate some companion 
flux in the process of attenuating speckles. To further calibrate our photometry, we introduce 
and measure the flux for fake point sources at random angles at a range of separations (0.5"- 
2.5") in each registered frame, rerun our ADI and LOCI pipelines, compute the attenuated 
flux in the final ADI/LOCI-processed images, and correct for this attenuation (~ 30% at 
GJ 758B's position for the LOCI reduction, ~ 10% for ADI). GJ 758B has L' = 15.97 ± 
0.19, where our photometric errors accounts for 1) the intrinsic photometric uncertainty 
for GJ 758B, 2) the intrinsic photometric uncertainty for our calibration star, and 3) the 
differences in photometry between the ADI and LOCI-reduced image^. Using the T09 H 
band photometry we then find H-L'= 3.29 ± 0.25 for GJ 758B. 

The H-L' color provides a first-order estimate for GJ 758B's temperature and spectral 
type (Table 1). GJ 758B is far redder than than expected if it were a L-T8 dwarfl. However, 
GJ 758B is comparable in color to Wolf 940B (H- L' = 3.38), a 3.5-6 Gyr- old, 570 K T8.5 
dwarf companion to a higher- mass brown dwarf (jBurningham et aL 2008). To derive GJ 
758B's effective temperature, we com pare its colors to th e colors predicted by the COND 
evolutionary models for its a ge range (iBaraffe et al.l 120031). Based on isochronal fi tting and 



age-activity correlations (e.g. iTakeda et al.l 120071 : iMamajek and Hillenbrandll2008l ). GJ 758 
has an age between 0.7 Gyr (~ 1 Gyr) and 8.7 Gyr (see d iscussion in T09). The 0.7 Gyr 



lower limit is based on isochrone fitting from the tracks ( iTakeda et al.l 120071 ). which may 
be highly uncertain given that GJ 758 is on/near the main sequence and given systematic 
disagreements between different sets of isochrones. For simplicity, we assume an age range of 
1 Gyr to 10 Gyr to cover the low end of the age range and the main sequence lifetime of a Sun- 



'"'If if our second, candidate point source is indeed the object T09 identify as "GJ 758C", its magnitude 
would be L' — 16.01 ± 0.38; the 3cr hmit at its location reported in T09 is ~ 16.03. Regardless of whether 
our candidate detection of "GJ 758C" is real or is a speckle, the H-L' color for the purported "GJ 758C" 
object must be bluer than w 2.5. 

^See the sample compiled by S. Leggett: [T]HREF]http://staff.gemini.edu/~sleggett/LTdata.html 
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like star (jSackmann et al.lll993l ). For this age range, GJ 758B has an allowed temperature 
range of log(Te) ~ 2.675-2.85 or Te = 560-9o K, where the extrema are determined from 
the ± la values for the H-L' color (3.04 and 3.54) and a slight age dependent calibration 
between temperature and color. 

The implied luminosity of GJ 758B as probed by the companion's absolute H band 
magnitude and H-L' color is broadly consistent with COND model predictions of GJ 758B's 
age range. For GJ 758B's H-L' color, the COND models at 5 and 10 Gyr predict an absolute 
H band magnitude of ~ 18.4 and 18.2, respectively. At 1 Gyr, the COND models predict 
Mh ~ 18.6, though given the uncertainty in the H-L' color the predicted near-IR absolute 
magnitudes are still consistent with observations. Given our range in effective temperatures 
and comparing L/Lq and Te from 1-10 Gyr for the COND models, the bolometric luminosity 
of GJ 758B is log(L/L0) ^ -6.I-0.2. 

Figure [2] plots the predicted temperature evolution of 10.5-42 Mj companions compared 
to the temperatures of GJ 758B and other cool, substellar-mass objects above and below 
the deuterium-burning limit. GJ 758B has an inferred mass ranging from 10-20 Mj (at 
1 Gyr) to ~ 25-40 Mj (at 10 Gyr), consistent with estimated mass for Wolf 940B (24-45 
M j). Within the context of the temperature evolution of substellar objects predicted by the 
COND mo dels, GJ 758B is likely higher mass th an IRXS J162041-210524B fILafreniere et al. 
20081 . boioh and 2M 1207B Jchauvin et al.ll2004h . GJ 758B is consistent with b eing an older 
cooler analogue to mid-L brown dwarfs in the Pleiades (e.g. BRB29) studied by lBihain et al. 
( 120101 . and references there in) and the recent ly discovered brown dwarf companion to PZ Tel 
a 12 Myr-old 1.3 M© star jBiller et al.lboioh . 



For most of the GJ 758B age range, the COND models yield estimated masses above 
the deuterium burning limit nominally separat ing planets from brown dwarfs. The more reli 
able age indicators, age-activity relations (e.g. iMamajek and Hillenbrandll2008l : iBarnes et al. 
20071 ). suggest that GJ 758 at least as old as the Sun, and perhaps as old as ~ 8.7 Gyr (T09). 
Thus, the mass of GJ 758B is most plausibly comparable to that of low-mass brown dwarfs, 
not high- mass planets. 

We investigate the range of plausible orbital parameters with a Monte Carlo simulation 
following the method of T09. With three astrometric points spanning one year (Figure [3l 
left panel), the orbital properties of GJ 758B are better constrained: only 1.2% of the trial 
orbital solutions "fit" the data compared to ~ 6% based on only the 2009 data (T09). The 
best-estimated eccentricity is high {eM.w. ~ 0.73~o-2i), and the best-estimated inclination 

+ 14° 

indicates that the orbit is viewed neither pole on nor edge on^iM.w. ~ 50°-22°). The best- 

+ 32AC7 

estimated semimajor axis is ~ 44 AU-i4ac7^ or 60% larger than its projected separation of 
~ 28 AU. 
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4. Discussion 

Our study recovers the detection of GJ 758B reported by T09 and demonstrates that it 
has a red H-L' color consistent with being an ultra-cool (T f» 560 K) low-luminosity com- 
panion to a solar-type star. This estimate is consistent with the low end of the temperature 
range (550-640 K) quoted by T09 based on assuming that the COND models accurately 
reproduce the fluxes of substellar mass objects in each near-IR bandpass, although the al- 
lowed temperature range is about the same as that quoted by T09, given our photometric 
errors. We also identify a second point source at a separation comparable to the putative 
"GJ 758C" companion reported by T09, though its detection has a low significance and it 
is likely a residual speckle. However, if this is the same object as that previously reported, 
it is not in a bound orbit around GJ 758. 

Based on the COND evolutionary models, GJ 758B has a mass between ~ 10 Mj and 40 
Mj. However, for most of GJ 758's age range, especially for ages derived from more reliable 
methods, its inferred mass is greater than the deuterium burning limit. Thus, based on our 
analysis, it is then much more likely that GJ 758B is a low-mass brown dwarf companion 
than a planet if "planet" is defined by an object's mass relative to the deuterium burning 
limit. 

A separate, perhaps more physically motivated criterion for defining an object to be a 
"planet" is its formation history: whether or not it formed in a proto planetary disk around 



a you ng star. Sophisticated numerical models for planet formation (IKenyon and Bromley 



20091 ) necessary for testing formation theories have yet to be directly applied to systems with 
directly-imaged planets. However, analytical arguments show that forming GJ 758B by core 
accretion is po ssible only if gas accretion is exceptionally efficient. Using Equation 14 from 



Rafikovl (120101 ). the minimum surface density of solids required to trigger core instability 



within th e lifetime of the dis k is ~ 6x10 ^ g cm ^ for a (long) protoplaneta ry disk l ifetime 



of 5 Myr (jCurrie et al.l 120091 ). This surface density is comparable to that for iRafikovi 's Min- 
imum Mass Solar Nebula profile at our best fit semimajor axis of 44 AU. Once a core has 
formed, the isolation mass to which a non-migrating gas giant could grow is > 10-20 Mj in 
disks containing that quantity of gas. However, accreting > 10-20 Mj of gas in a disk that 
cannot be more than ~ 100 Mj in mass total (else the disk would be gravitationally unsta- 
ble) requires the planet to accrete very efficiently. Companion formation by gravitational 
instability of the disk, in contrast, can only occur while the disk is still accumulating gas 
from the protostellar core, and preferentially for ms companions with masses rnore s imilar to 
field brown dwarfs than jovian-mass planets fe.g. iKratter et al.lboiol : iRafikovlboOsI l 



Formation may also be possible earlier, during the fragmentation of the molecular cloud 
core that formed GJ 758: GJ 758B would then comprise the low-mass end of the binary 
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mass function. Core fr agmentation can pro bably produce objects much less massive than 
GJ 758B (e.g. 5-10 Mr lTodorov et al.ll201Cll ). so we consider this scenario to be a plausible 
one for forming GJ 758B in addition to fragmentation during the protostar phase. Future 
multiwavelength studies of GJ 758B will better constrain its temperature and atmospheric 
properties. Combined with improved age estimates derived from stellar rotation and ac- 
tivity, these data will better place GJ 758B within the context of directly imaged low-mass 
brown dwarfs useful for investigating the mass function of objects approaching the deuterium 
burning limit. 
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Fig. 1. — LOCI-reduced images from our May dataset (left) and June dataset (right) shown 
at high contrast to display residual speckle noise. GJ 758B (lower circle) is recovered in 
both data. A candidate point source with a separation comparable to that reported for "GJ 
758C" (top circle, left panel) is also detected in the more sensitive May data, though this 
may instead be a residual speckle. Because of shorter integration times and far poorer field 
rotation, residual speckle noise is far more severe for the June data, and we do not recover 
the second, candidate point source. 
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Fig. 2. — Effective temperature vs. age for low-mass brown dwarfs/planetary- mass 
objects orbiting solar and subsolar-mass stars. Along with GJ 7 5 8B, these include 



Chauvin et al. 



PZ Tel teiller et all 120101) . IRXS J160929.1B flLafreniere et al. 



20041). t wo L dwarfs in the Plei ades (iBihain et al. 



flBurningham et al.ll2008h . and SDSS J1416-I-13AB flBurningham et al. 



2008 



20101), 2M 1207B 
20ld\ . Wolf 940B 
20091 ). We overplot 



the temperature ev olution for 10.5-42 Mj 
jBaraffe et al.lboosi ). 



objects from the COND evolutionary models 
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Fig. 3. — (left) Projected separation for GJ 758B between May 2009 and May 2010, com- 
bining our data (upper right) with that from T09. The arrows show the motion of GJ 758B 
from May 2009 (lower left) to May/ June 2010 (upper right); the diamonds are the best- 
fit positions based on a rectilinear model, (right) Range of orbital solutions for GJ 758B 
following the method of T09. 



- 14 - 



Table 1. Observed and Derived Properties for GJ 758B 



Pro]:)erties 


GJ 758B 


Photometry 




App. V 


15.97 ± 0.19 mag 


Abs. L' 


14.84 ± 0.21 mag 


H-L' 


3.29 ± 0.25 mag 


Derived Physical Properties 




Approx. Spectral Type 


T8.5-T9 


Te (inferred) 


^ 560 +150, -90 K 


log(L/L0) (inferred) 


-6.1 +0.3,-0.2 


Inferred Mass (Mj) (1 Gyr) 


10-20 


(5 Gyr) 


21.0-31.5 


(10 Gyr) 


25-40 


Astrometry 




May 3, 2009 (T09) 




Separation (") 


1.879 ± 0.005 


Position Angle(°) 


197.17 ± 0.15 


August 6, 2009 (T09) 




Separation (") 


1.858 ± 0.005 


Position Angle(°) 


198.18 ± 0.15 


May 27, 2010 (this work) 




Separation (") 


1.823 ± 0.015 


Position Angle(°) 


199.76 ± 0.15 


Projected Separation (AU) 


28.36 


(5(Sep., ") 


0.035 


J(PA, <>) 


1.58 


vr 


-0.9 ± 0.2 AU/yr 




0.94 ± 0.1 AU/yr 


Derived Orbital Properties 




Semimajor Axis (AU) 


44.12 (30.28, 76.66) 


Eccentricity, Inclination 


0.73 (0.52, 0.85), 49.64° (27.87°,63.34°) 



Note. — The changes in separation (5Sep.) and position angle ((5PA) are given between 
our May 27, 2010 data and the August 6, 2009 data presented by T09. The semimajor axis, 
eccentricity, and inclination listed refers to the median-weighted value and the two values 
defining the inner 68% of weighted values. The instantaneous radial and angular velocities 
(v_R, ^pa) were determined at an epoch of 31 Dec 2009 (the error-weighted mean of our May 
datum and T09's data) and an assumed distance of 15.5 pc. 



